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Abstract 

Certain quantum topological invariants of three manifolds can be written in the form 
of the Gaussian sum. It is shown that such topological invariants can be approximated 
efficiently by a quantum computer. The invariants discussed here are obtained as a 
partition function of the gauge theory on three manifolds with various gauge groups. 
Our algorithms are applicable to Abelian and finite gauge groups and to some classes of 
non-Abelian gauge groups. These invariants can be directly estimated by the nuclear 
magnetic resonance (NMR) technique used for evaluating the Gaussian sum. 
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1 Introduction and Summary 



Despite of the expectation for potentially wide applications of quantum computation[l], the 
speedup due to quantum algorithms is limited only to restricted problems. Well-known ex- 
amples with exponential improvement such as the factorization and discrete-log algorithm[2] 
both depend on the quantum Fourier transformation as a crucial technique. The similar 
method was later extended to other number theoretic algorithms[3, 4]. 

Experimental implementations of these algorithms are even more limited due to the dif- 
ficulty of maintaining coherence during quantum computation even for a small number of 
qubits. Among them, the most remarkable example is a realization of the Shor's factorization 
algorithm [5] by the NMR technique. The factorization of larger numbers was also demon- 
strated recently by the NMR technique [6] through evaluating the Gaussian sum. These 
methods crucially rely on the experimental determination of periodicity in a function. It is 
worth investigating to find more applications based on similar techniques. 

In this work, we explore different applications of the Gaussian sum through its relation 
to the topological invariants of three dimensional manifolds. The evaluation of topologi- 
cal invariants of three manifolds has been discussed in relations with topological quantum 
computation [7, 8]. A similar line of thoughts lead to efficient algorithms for the estima- 
tion of topological invariants. The efficient algorithm approximating Jones polynomial was 
constructed [9, 10, 11]. Although these approaches are theoretically elegant and intriguing, 
experimental realization of these algorithms is still far from a practical stage. 

Quantum topological invariants can be obtained essentially as partition functions of some 
special classes of quantum fields defined on three manifolds [12]. The further simplification 
may be possible to reduce the partition function to the form with its periodicity manifest. 
In some cases, it can be reduced to the Gaussian sum form. 

In this paper, we discuss three cases in which this reduction is possible. Each case is 
related to an interesting physical theory. In Sec. 2.1, we discuss Abelian gauge theory 
invariants. They are defined as partition functions of the Abelian Chern-Simons gauge field 
on the underlying three manifold. A large class of such invariants can be written as Gaussian 
sums. In Sec. 2.2, we discuss non-Abelian gauge theory invariants, defined as partition 
functions of the non-Abelian Chern-Simons gauge field on three manifolds. We study the 
case in which these variants can be reduced to Gaussian sums.. Finally in Sec. 2.3, we discuss 
the case in which the gauge group is finite. In this case, we obtain the invariants known 
as Dijkgraaf-Witten invariants[13]. In these cases, we show that efficient approximation 
algorithms based on the evaluation of multivariate Gaussian sums exist. These invariants 
can be experimentally observed by the NMR technique for factorizing large numbers and by 
other optical processes. 
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2 Wilson loop as a topological invariant 



The Chern-Simons action for the SU (2) gauge field A is 

k 



AidjAk + ^ AiAjAk 



(1) 



S[A] [ (fxe'^'^Tr 

An Jm 

The partition function we study is 

Z[M] = [ DAe'^^^^l (2) 

JM 

The observable we study is the Wilson loop {W{C)) defined along an oriented curve C as 

W{C) =TTR[Pe''fc^^'^^'], (3) 

where the trace is taken in the irreducible representation of a gauge group. By closing a loop 
and taking a trace, W{C) is defined to be gauge-invariant. When the curve C is a hnk L 
consisting of disjoint knot components Ci, C^, 

m 

W{L) = WiC^, ...,C^) = l[ WiC,). (4) 

The expectation value of the Wilson loop is given by 

{W{L)) = [ DA W{L)e'^^^\ (5) 

JM 

Since {W{L)) docs not involve any metric, we expect that it is essentially topologically 
invariant. In order to see the exact topological invariance, we consider the Dehn surgery 
representation of the three manifold M with a framed hnk L embedded in as follows: 
First consider a torus obtained as a tubular neighborhood of each component of L. We 
denote a union of all m tubular neighborhoods as N(L). Wc can decompose as a union 
between S*"^— int N{L) and m tori. The Dehn surgery is defined by pasting each torus back 
to where it was in with a twist such that a meridian of each torus and each longitude of 
5"^— int N{L) is identified. Arbitrary closed oriented three manifold M can be obtained by 
Dehn surgeries[14] around a hnk L embedded in S^. Thus regarding a twist as a framing 
of each component of a link, topological information of M can be obtained by studying a 
framed link in S^. 

Two three-manifolds obtained by different framed links are homeomorphic to each other 
if and only if they are related by Kirby moves[15]. In particular, any topological invariant 
needs to to be invariant under Kirby moves. Kirby moves consist of adding or removing a 
trivial knot C with its self- linking number ±1 and making a connected sum of a component 
of a link with another along their framing. It turns out that the expectation value {W{L)) 
is not invariant under Kirby moves but aquires +1, a phase as 

{W{L))^e^^'{W{L)), (6) 

where 9 = 3n{k — 2) /4k (assuming A; > 2) is related to the central charge in two dimensional 
theory. Thus {W{L)) itself is not a topological invariant and needs some modification. We 
will come back to this in Sec. 2.2. 
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2.1 Abelian gauge group 

We construct the topological invariant for the Abehan gauge field A in the following. For 
the Abelian gauge field A defined on the three manifold M, the second term in (1) vanishes 
and its action S becomes 



S[A] = ^ [ d'xe'^^AdjAk. 



(7) 



We consider a hnk composed of m components Ci, ...,Cm embedded in M = 5'^. The 
vacuum expectation value of the Wilson loop is related to the correlation function for the 
gauge field as 



{W{C^,...,Crn)) 



i=i 



= 1 



= exp 



m r r 

^iCj j> j> {Af,{xi)A,,{xj))dxfdxj + ■ 

ij=l JCj 

m r r 

^ dCj j> j> {Af,{xi)A^{xj))dxfdx'< 

• 1 'J Ci " C-i 



The correlation function of the Wilson loop can be written as 



{A^{xj)A,{xi)) = ^ / dx^ [ dyh 



\x — 



(8) 



(9) 



If we write the chare; integer multiples of an elementary charge e as Cj = njC (n^ = 

1,2,...) and, by absorbing e in the coupling constant as k = ko/e'^, then the correlation 
function has the following form [17]: 



where 



Jij — - 



— / dx' [ dy^t'^^ 
2n JCi JCi 



\x — 



(10) 



is a linking number of links Cj and Cj and is a self-linking number. Note that is a 
symmetric integer matrix and is an integer vector. Naive definition of 



(12) 



4 



contains divergence and we need to regularize it. In order not to sacrifice the metric inde- 
pendence by explicitly introducing the cutoff parameter, we choose the point splitting regu- 
larization prescription by using framing of links, i.e. deforming the contour C parametrized 
as x^{s) to the new contour C5 by x^{s) — > x^{s) + 5x^{s) and define 



lim — I dx' I dyh'^'' ^^ (13) 
Sx^o2nJc Jc, \x-y\^ ^ ' 



which is a self-linking term between the original knot and its framing. 

The Wilson loop is invariant under the shift of variables — > + /c for an integer k. 
Hereafter we only consider this case. This will limit the range of rii to be finite: rii = 1,2,..., k, 
where each rii gives a different representation of f/(l). In order to obtain the quantity 
independent of the representation, we sum over all representations and write this sum as 
{W{L)) as follows 



(14) 



ni,nj=l 

If we write Ju — rsi in above. 



(15) 



This can be viewed as a partition function for the spin variables rii that take k different 
values at the vertices i of the graph. The shape of the graph is determined by the edges 
connecting vertices i and j for Jij 7^ 0. 

We assume that k = p™' with an odd prime p, Jij can be diagonalized modulo k with a 
matrix U e SL{m, Z)[18]. Let us denote the components after diagonalization as Ji, Jm- 
Then U'^JU = ©*Li Ji and 



i=l \ni / 

For k = 1 (mod 4), {W{L)) itself can be shown to be a topological invariant [19]. In this 
case, the Abelian invariant ta{M) can be written as a product of Gauss sums: 



TA{M) = llG{k,Ji), (17) 

where 



1=1 



k-l 



G{k,a) = J2e~'^ ■ (18) 

n=0 
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Appearance of G{k, a) indicates that ta{M) is computationally hard to evaluate as the 
estimation of the Gaussian sum is considered to be classically hard. On the other hand, the 
use of quantum algorithms enables us to approximate ta(M) in polynomial time, as each 
Gaussian sum factor can be done so[4]. 

While the norm of the sum is straightforward to calculate, its phase requires 0(l/e) times 
measurements to estimate within the error e. The Gauss sum is defined by 

G(A;,a) = X:xWe-'^"" (19) 

n=0 

x(n) is a multiplicative character. For our purpose, we take x(n) to be the Legendre symbol: 

X{n) = (^) . (20) 

It is convenient to define x(0) = 0. For an odd prime k, 

= n('=-^)/'(mod k) for n ^ (21) 

holds. One can easily show that G{k, a) for an integer a is equivalent to the familiar form 
given in (18) [20]. 

The Gauss sum can be efficiently estimated by the following steps[4]. First we construct 
a quantum state with the coefficient given by the Legendre symbol. Second we use a unitary 
transformation to this state to change the coefficient to the Gauss sum. Then we make a 
measurement to obtain the amplitude of the state closely related to this state that gives the 
estimation of the Gauss sum. 

Preparation of a state: 

\x)^^^j:xin)\n) (22) 



can be done in the following way. We use the fact that the following quantum Fourier 
transform 

b)^^X:e-xP^|.) (23) 

can be performed efficiently on a quantum computer by 0(A;^) steps[l]. Making the quantum 
Fourier transformation on the second term on the product state |n)|l) gives 

-I k — l 

N^^T^E^-IO, (24) 
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where q = k . By shifting each \l) by ^ log^ n 

^~^\og.n) (25) 



= <1 



\n) 




fc-1 ( 


\n) 




fc— 1 


^ logq n 


\n) (g) 


1 fc-1 




\n) <S) 


-1 fc-1 



where the base q of the logarithm is viewed as the element in Z^. log^n is an integer p 
such that — n(modA;). The evaluation of the discrete logarithm log^n can be efficiently 
performed on a quantum computer by poly(log2 k) steps[2]. 

We make another quantum Fourier transform and write the right hand side in (22) as 



1 1 A;— 1 fc— 1 

Vk-lVk „=o 
1 G{k,a) '^ 

Vk^ Vk ti 



(26) 



where we used 



r(h nn - X{1 ^)G{k,a) for/^0 . . 

G[k,al)-^ for/ = ^"^^^ 

from the first line to the second line. Mapping \l) to x(/)^IO by following the same procedure 
as in (26) will make the last term in (26) into 

^|X>. (28) 

Knowing that \G{k,a) \ = Vk[20], now the phase information for the Gauss sum is stored in 
the coefficient. 

By operating a phase gate which adds the phase G{k, a)/\/k — e"^'^ on the state |x) by 
the above procedure and does nothing on the state |0), we can transform the initial state 
|0) + Ix) to 

|0) + |X) - |0)+e-^|x). (29) 
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Then by making observation in the order of 1/e times, we can determine the phase within 
the error e. 

Using the algorithm explained above, evaluating each Gaussian sum appeared in our 
topological invariant 

m 

TAiM) = l[Gik,J,) (30) 

i=l 

takes poly (fc, 1/e) time. Thus the evaluation of the invariant can be done in m x poly (A;, 1/e) 
= poly(m, fc, 1/e) time. 

2.2 Non-Abelian gauge group 

{W{L)) can be considered naturally as a link polynomial. For M = S*^, it is possible to 
choose the framing such that all self-linking terms vanish (called standard framing). In this 
case, {W{L)) only depends on the link L and becomes proportional to the Jones polynomial 
J(L, q). {W{L)) can be calculated explicitly by the combinatorial method similar to the one 
for the Jones polynomial. 

X 

LO 

Figure 1: The skein relation 

By assigning the same color n = 2 (corresponding to two-dimensional irreducible repre- 
sentation of sl2{C)) for all link components, above {W{L)) with local deformations shown 
in Fig. 1 obeys the following skein relation 

q'/\W{L^)) - q-'/\WiL^)) = {q^'^ - q-^"){W{L,)) . (31) 

W{L) still depends on the framing. In order to cancel frame dependence, we define {W{L)) 
by multiplying {W{L)) with [2]~-'^e~^'^™*^^^/^*', where [n] = sm{mi/k) / sin(7r/fc) and w{L) is 
a writhe defined by a difference of the number of positive and negative crossings. Then we 
see that {W{L)) satisfies 

q{WiL^)) - q-'{W{L_)) = (g^^ - q-'/'){WiLo)) . (32) 

Comparing this with the skein relation satisfied by the Jones polynomial J{L) 

t-^J{L+) - tJ{L^) = (t^/2 - r^/') J(Lo), (33) 
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we sec {W{L)) ioT n — 2 is essentially the Jones polynomial J{L+) by identifying t^/^ = 
— Thus we define J{L), a version of a Jones polynomial proportional to {W{L)) by 
normalizing J{L) as J{U) — 1 for a trivial knot U. The explicit relation between these is 
obtained in [21] as 

{W{L)) = [2]e^'^'/''=^r.=i-^*^J(L). (34) 

Mathematically rigorous definition of the invariants is given in [22]. It can be expressed 
in our formula (??) as follows. For each link component Q whose framing r^, we sum over 
all the color rii with proper weights. By writing the color for each link and the sum explicitly 
in (??), we have 

TNAiM) = C^e-^'''^ Y.[n^]-[nm]{Wn,...nM). (35) 

where c = (2/A;)^/^sin (tt/A;). Here we followed the definition in [21], which is different from 
that in [22] by an overall factor c'^, where v is the first Betti number of M. 

The evaluating tna{M) for general k is difficult. Neverthess, for certain values of k, the 
calculation can be simplified. In particular, for A; = 3, [2] = 1 and J(L) = 1 in (34) gives 

In [21], this result is used to write the invariant as 

tna{M) = 2-™/2g-*-aL/4 ^ g.-E,,,es-^i^/2^ (37) 

SdL 

where the sum is over all the sublink S with color n = 2. It is easy to see that this can also 
be written as follows: 

2 

tna{M) = 2-'^/^e-''"'^/^ ^ (38) 

ni=l 

Thus it is reduced to the multivariate Gaussian sum. From our result in Sec. 2.1, this sum 

can be estimated in m x poly(l/e). Combining with the calculation 0"^, that takes typically 
O(m^) time, we see that TjqA{M) can be estimated in poly(m, 1/e) time. 

2.3 Finite gauge group 

Geometrically the gauge field can be viewed as a connection on a principal G bundle over 
M. Then the path integral is a sum over all connections on M x G with the gauge group G. 
For a compact simply-connected gauge group, a principal G bundle is topologically trivial 
and the path integral over connections is reduced to the ordinary path integral of the gauge 
field. In our case, the action is given by the Chern-Simons action in (1). 
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For a finite gauge group G, a principal G bundle has a unique flat connection and nontriv- 
ial contribution to the path integral is coming only from a sum over the paths with different 
topologies. In other words, the path integral is replaced by the discrete sum over different 
conjugate classes of the G-bundle. Each conjugate class of the G-bundle is represented by the 
holonomy 7 for the map ni{M) — > G. In [13], the partition function for the gauge field with 
a finite gauge group is considered. The path integral is given by the sum over all homotopy 
classes of 7 as 

Tw(M) = ^^e2-^^ (39) 
l<^l 7 

where e'^'''^'' = (7*0, [M]) is a pairing between the three-dimensional homology class of M, 
[M] e H2.{M,U{1)) and the cohomology class 7*0; for each a e H^{BG;U{1)) pulled back 
by regarding 7 as a map from M to the classifying space BG. For a flat connection, the 
principle G bundle is completely specified by the homotopy class of 7. 

For G = Zk, more explicit expression can be obtained[18]. For the action S*^ to be 
real, {'j*a,[M]) G U{1). Since H^{BZk\U{l)) = Z^., choosing a is equivalent to choos- 
ing an integer / G 1, ...,k and e^'^*'/'^' G U{1). Meanwhile 7 has a corresponding element 
70 e H\M-Zk) by using Hom(7ri(A/), Z,,) ~Hom(i/i(M; Z), Z^) ~ H\M;Zk). Let 5* 
be a connecting homomorphism between H^{M;Zk) to H'^{M;Z) associated with an ex- 
act sequence O^Z^Z^Zk^O. Then 70 U 5*7o defines an element in H^(M; Zk). 
Exponentiating this will give an element in H^{M,U{1)). Coupling this element with [M] 
gives an action e^""''^^' = (7*0, [M]) = (e2'r*'(7ou5*7o)/fc^ [jvf]). 70 U5*7o can be evaluated as the 
intersection number among dual chains and given as = J27^j=i Jij'ninj/k[18], where the 
linking number Jij appears when M is represented by Dehn surgeries along the link. Thus 
we obtain 

TDwiM) = ^ E e'^^^^r.- '^--'-^■/^ (40) 

where we took / = k. Thus again the invariant TDwiM) can be written as the multivariate 
Gaussian sum, thus can be estimated efficiently with the method developed in Sec. 2.1. 
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